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A field configuration utilizing local static and oscillating fields is constructed to achieve universal 
(but low-order) protection of two-qubit states. That is, two-qubit states can be protected against 
arbitrary system-environment coupling with a driving field whose frequency is sufficiently large as 
compared with the cutoff frequency of the environment. Equally important, we show that it is 
possible to construct driving fields to protect two-qubit entangling gates against decoherence, with- 
out assuming any particular form of system-environment coupling. Using a non-Markovian master 
equation, we further demonstrate the effectiveness of our continuous dynamical decoupling fields in 
protecting entanglement and the excellent performance of protected two-qubit gates in generating 
entanglement. The results are complementary to current studies of entanglement protection using 
universal dynamical decoupling pulse sequences. 



I. INTRODUCTION 

The unwanted interaction between a system and its en- 
vironment causes decoherence, i.e., the loss of quantum 
coherence. Since future quantum technologies rely on co- 
herent quantum states, it has become increasingly impor- 
tant to effectively suppress decoherence. To achieve this 
goal, various schemes, such as error-correction codes 
decoherence-free subspaces 0, ^ and dynamical decou- 
pling (DD) have been proposed. It is expected that 
in a large-scale working quantum computer, all of these 
schemes will be used in some way to store quantum states 
with high fidelity. For instance, a combination of dynam- 
ical decoupling and quantum error-correction codes has 
been proposed to combat errors due to spontaneous emis- 
sion (see also Ref. for a recent study on such a 
hybrid scheme). 

Our focus here is on DD. In the DD approach, exter- 
nal time-dependent fields are applied to the system such 
that the interaction term between the system and its en- 
vironment rapidly flips sign. In this way, the effect of the 
environment on the system is canceled to a certain de- 
gree. The key advantage of DD, compared to some other 
methods such as quantum error-correction codes, is that 
no overhead is required - the qubits storing quantum in- 
formation are protected directly, without any need for 
extra qubits. Moreover, DD requires neither quantum 
measurements nor feedback control. 

Broadly speaking, two types of DD have been studied: 
pulsed DD, which uses sharp pulses (impulsive pulses in 
many cases) to counter the effect of the environment, and 
continuous DD. which uses continuous- wave driving, that 
is, fields with simple harmonic time dependence. We dis- 
cuss pulsed DD first. The pioneering work on pulsed DD 
considered pulses applied, at equal time intervals, to a 
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single qubit coupled to an environment Many ex- 
tensions have been worked out since then, including, for 
example, DD for an arbitrary finite-dimensional system 
Q and the suppression of arbitrary internal coupling in a 
quantum register [llj . Studies of non- impulsive DD are 
also fruitful, by use of, for example, the so-called "Eule- 
rian DD" schemes [ij, fl3| . a nd optimized pulses under 
an energy cost constraint |I4{ or a minimum leakage re- 
quirement [iHl • Recently, a significant advance was made 
when it was shown by Uhrig [16| that by using aperiodic 
pulses, the so-called Uhrig's DD (UDD) scheme, the co- 
herence of a single qubit can be protected to the iV*^ 
order by using only N (or iV -|- 1) instantaneous pulses. 
Uhrig originally considered only pure dephasing in the 
spin-boson model. Yang and Liu |I7| then showed that 
UDD is universal in the sense that it does not depend 
on how a single qubit is coupled to its environment. Go- 
ing further. West el al. [l^] constructed a nested UDD 
sequence that can protect a single qubit against both de- 
phasing and relaxation at the same time. A mathemat- 
ical proof for the effectiveness of nested UDD sequences 
has been recently given in Ref. [l^. Effects of nonideal 
pulses on UDD are also under investigation [2^ . On the 
experimental side, the excellent performance of UDD in 
protecti ng sing le-qubit quantum states has been studied 
in Refs. [SHII- 

Given high-efficiency single-qubit DD schemes, extend- 
ing single-qubit DD to two-qubit (or multi-qubit) deco- 
herence control becomes more interesting. It should be 
emphasized at this point that two-qubit (multi-qubit) de- 
coherence control offers a whole new set of challenges (see 
also 0, [ill[i3)- For instance, a fundamental objective 
of two-qubit decoherence control must be to protect two- 
qubit entanglement, since entanglement has been iden- 
tified as the key resource for quantum information [2^ . 
As is now well known, quantum entanglement can, unlike 
single-q ubit coherence, vanish in a finite amount of time 
25l-l30l| . In addition, there are different types of system- 
environment coupling that are not present in single-qubit 
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cases. In particular, there can be a noisy interaction be- 
tween two qubits, and errors such as correlated bit flip- 
ping, dissipation, and dephasing might emerge. Along 
this general direction of two-qubit DD, Ref. ^Jj showed 
for the first time that it is possible to construct a pulse 
sequence to protect a known two-qubit quantum state to 
the A'^th order using N pulses, without any knowledge 
of system-environment coupling. During the same year, 
the same group of authors advocated the use of a nested 
sequence of UDD pulses to protect unknown two-qubit 
states with high efficiency, with each layer eliminating 
different noise terms |32| . It is now clear that to protect a 
completely unknown two-qubit state with high efficiency, 
four layers of UDD pulses are required [sJ^lH]. These 
schemes are exciting because they are universal. That is, 
so long as the pulses arc applied fast enough (as com- 
pared with the cutoff frequency of the environment), we 
do not need to assume anything about the actual form of 
the system-environment coupling. Nested-UDD schemes 
have also been extended to multi-qubit systems with re- 
markable mathematical insights [H, [sj] . Parallel with 
these theoretical advances, preliminary experiments on 
entanglement protection using pulsed DD have been per- 
formed in Rcfs. [31433. 



Interestingly, many pulsed DD schemes mentioned 
above can be considered to be too strong in the sense 
that, while it does protect a quantum state with high ef- 
ficiency, it also generally freezes useful coherent evolution 
generated by the system's own Hamiltonian. To achieve 
useful coherent evolution concurrently with pulsed DD 
[H, [11] , one idea is to encode the logical qubits in phys- 
ical qubits and then design the fields in such a way 
that the gate operation commutes with the pulse opera- 
tions However, apparently an overhead is required. 
Recently, a general procedure utilizing finite-power and 
finite-bandwidth pulses has been worked out for con- 
structing dynamically corrected gatesfDCG) without en- 
coding or measurement overhead [4l|, [4^. Going fur- 
ther, by concatenating DCG's, it is possible to achieve 
arbitrary accuracy in quantum gate implementation [43j . 
However, for arbitrary system-environment coupling, dy- 
namically corrected two-qubit gates have not been explic- 
itly constructed and it is unclear how complicated the 
solution might be. 

The existence of universal UDD schemes for two-qubit 
entanglement protection motivated this work. In par- 
ticular, we shall investigate the usefulness of universal 
continuous DD (more specifically, DD based on driving 
fields with simple harmonic time dependence) to protect 
two-qubit states, and hopefully, also to protect two-qubit 
gates. Switching from pulsed DD to continuous DD, the 
sacrifice is obvious as compared with UDD and nested- 
UDD schemes: the performance of continuous DD is of 
a low-order nature. But our interest here is not with 
the high-order performance of a DD scheme. Rather, we 
ask the following important question: Are there univer- 
sal continuous DD schemes to protect two-qubit states 
and two-qubit gates, irrespective of how a two-qubit sys- 



tem is coupled with its environment? This is a perti- 
nent question to ask because, compared with pulsed DD, 
continuous DD has some advantages from a practical 
point of view. For example, there is no longer any con- 
cern about pulse timings or pulse-sequence engineering, 
and the higher driving frequencies that wc can achieve 
with continuous fields arc naturally expected to elimi- 
nate higher frequency noise sources (43 |. Indeed, contin- 
uous DD schemes to protect a single qubit have attracted 
considerable interest j43 - [47| . Among the known features 
of single-qubit continuous DD, most relevant here is the 
fact that continuous control fields may be constructed 
to protect a quantum state and implement agate at the 
same time, without the use of any overhead [i^, [13] (thus 
forming a type of DCG [4l|, [12] )■ However, it is im- 
perative, considering the complexity of two-qubit deco- 
herencc, that such Hamiltonians be constructed for two- 
qubit gates as well. After all, in the circuit model of 
quantum computation, two-qubit gates are of fundamen- 
tal importance. 

What is lacking currently is a completely general treat- 
ment of continuous DD for two-qubit systems. A recent 
study [i^ considered the application of external fields to 
protect a multi-qubit system against a restricted class of 
dephasing and relaxation mechanisms. Some stimulat- 
ing progress has also been made in Ref. [11] . However, 
therein only local noise terms were considered, which 
amounts to making a specific assumption applicable to 
only one class of system-environment coupling. With 
such an assumption, continuous fields were constructed 
for the protection of two-qubit quantum states and two- 
qubit gates against decoherence [i^. Nevertheless, as 
also seen below, if only local noise terms are considered, 
then the issue of decoherence control is somewhat quite 
analogous to singlc-qubit continuous DD and is not uni- 
versal. The explicit task of this work is hence to ex- 
tend the work in Ref. [i^ to cases with the most gen- 
eral system-environment coupling. Our universal two- 
qubit continuous DD schemes presented below may be of 
great use if very high control fidelities are non-essential. 
Certainly, universal continuous DD schemes may be also 
combined with pulsed DD for hybrid DD schemes. 

The organization of this paper is as follows. In Sec. II, 
we construct continuous driving fields that are able to 
protect an arbitrary two-qubit quantum state against its 
environment. Once we have such a scheme, we use this 
knowledge for the construction of two-qubit gates pro- 
tected against the most general decoherence. It is shown 
in Sec. HI that an explicitly constructed control Hamil- 
tonian can implement and protect two-qubit gates at the 
same time, and can therefore yield much better gate per- 
formance. In Sec. IV, in order to test these theoretical 
results, we introduce a non-Markovian master equation. 
By modeling the environment as one (or many) thermal 
bath(s) possessing an Ohmic spectral density, we show, 
in Sec. V, the results of our numerical simulations illus- 
trating the excellent performance of our fields in protect- 
ing two-qubit states as well as in implementing two-qubit 
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gates. Finally, Sec. VI concludes this paper. 



where 



II. UNIVERSAL CONTINUOUS DYNAMICAL 
DECOUPLING 

We start off by considering the Hamiltonian of a two- 
qubit system interacting with its environment (modeled 
by - possibly more than one - thermal bath later), 



Hq + Hb + HsB 



(1) 



where Hq denotes the Hamiltonian of the two-qubit sys- 
tem, Hb the Hamiltonian of the environment, and Hsb 
is the interaction Hamiltonian between the system and 
its environment. 

The most general form of Hsb is given by [lH 



H 



SB 



3 3 
3(1)^(1) 



k=l 



k=l 



- E 

k=i,i=i 



(12)^(1)^(2) 



(2) 

where ci = ax,cr2 = <Jy,a3 = cr^ , and the B opera- 
tors denote arbitrary environment operators. Note that 
this form is considerably more complex than the local- 
environment interaction Hamiltonian considered in Ref. 



since nonlocal terms such as B 



(12)^(1)^(2) 



kl 



taken into account. 

We now consider continuous driving fields applied to 
the system, whose effect is described by the Hamiltonian 
Hc{t). Corresponding to Hc(t), there is a unitary opera- 
tor Uc{t), given by the time-ordered exponential of Hc{t) 
(h ^ I and T is the time-ordering operator throughout), 



Uc{t) = r cxp 



—I I Hc{s) ds 





(3) 



In order to achieve continuous DD, Uc{t) must fulfill two 
criteria. The first is that it should be periodic in time 
with a period denoted by tc, that is. 



Uc{t + tc) = Uc{t). 



(4) 



Secondly, in order to decouple the system from the envi- 
ronment, we hope to have 



Ul{t)HsBUc{t) dt^O. 



(5) 



Technically, these conditions can be derived using the 
Magnus expansion. For completeness, following the 
treatment given in Ref. jsoj . we show that these condi- 
tions indeed lead to a low-order decoupling of the system 
from the environment. Ideas and notation introduced 
here will be used again when we explain the reasoning 
behind the construction of control fields for gate protec- 
tion. 

The Hamiltonian for the total system in the presence 
of the control fields can be written (in the 'lab' frame) as 

Htot = Ho + H,{t) +Hb + Hsb = H' + H,{t), (6) 



H' — Hq + Hb 



H 



SB- 



(7) 



Our goal is to see how a state evolves under the action 
of this total Hamiltonian, if conditions (|4]) and (O are 
satisfied. In order to do so, we transform to the frame of 
the control fields, that is, we rotate the basis by Uc{t). 
Then, in this frame, a total system-environment state 
evolves under the action of the unitary time-evolution 
operator. 



Utotit) = Texp 



-i / H'{s)ds 
Jo 



(8) 



where H'{s) = Ul{s)H'Uc{s). 

At time t ~ Ntc {N is a positive integer), because 
H'{s) is periodic with period tc, we have 



Utotit) = Utotitc) 



N 



and 



Utotitc) = Texp 



-I I H'is)ds 





(9) 



(10) 



The Magnus expansion [5l| allows us to write (refer to 
Appendix A), 



Utotitc 



exp 



-zt,(i/(") + (1) + . . .)' 



with, 



H(0) ^ 1 / ds H'is). 



(11) 



(12) 



We ignore the higher order terms since we are concerned 
with a low-order DD only. 

Now it is at this point that the condition Eq. ([5]) comes 
in. Because of this condition, i.e., J^" Ulit)HsBUcit)dt — 
0, we can eliminate the Hsb term in H^^\ Therefore, 

nm = h + Hb, with 



H^L r dsUlis)HQUcis). 

tc Jq 



(13) 



We also note that H is independent of tc- Since Udt) is 
periodic in time with period tc, we can write it as some 
function oft/tc, say, U'^t/tc)- Then, 



1 



H = -l dtU'Jit/tc)HQU'dt/tc), 
c Jq 
1 

dx U'Jix)HoU'dx), 





(14) 



with X = t/tc- It follows that H is indeed not an explicit 
function of tc- 
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Keeping in mind that N = t/tc, wc then find that to 
lowest order in i^i 



Utotit) 



t/ta 



(15) 



Finally, transforming back to the original frame (the 'lab' 
frame), we find that the unitary evolution operator in this 
frame is 



Utotit) w Uc{t)e 



(16) 



But, because of the condition t = Ntc and the period- 
icity of Uc{t), Uc{Ntc) is just identity. Obviously, then, 
the system has been decoupled from the environment - 
they both evolve independently, since H acts only on the 
system Hilbert space, while Hb acts only on the environ- 
ment Hilbert space. Roughly speaking, we can under- 
stand this result by realizing that under the condition in 
Eq. ([S]), the system-environment interaction is averaged 
out in the frame of the control fields. 



A. Suppression of local noise 

Let us now come back to our problem of finding con- 
trol fields to protect an arbitrary two-qubit state. As 
stated before, in Ref. [1^, continuous fields were found 
to eliminate the local noise terms. Local noise here means 
the following restricted form of system-environment cou- 
pling. 



Hi 



SB 



k=l 



(1)^(1) 
k '^k 



(2)^(2) 
k ■ 



(17) 



Other coupling terms in Eq. ^ different from above are 
loosely called nonlocal noise terms. 

We recap what is already known - how to find con- 
tinuous fields to eliminate the above-defined local noise 
terms, as is done in Ref. [i^ (but with more details). We 
first observe that with the unitary control operator 



where 



U.{t) = U^'\tM'\t), (18) 



U^^''\t) = e-2'^*'^i""-*/tc^ fc = 1,2, (19) 



with Tlx a non-zero integer, we eliminate noise terms pro- 
portional to CTy'^' and cri'^^. Intuitively, this follows from 
the fact that Uc{t) is just a rotation operator, and there- 
fore causes the Cy'"^ and cri'^'' noise terms to rotate so that 

(k) 

they average out to zero. However, it leaves the cri noise 
terms untouched. In order to cancel these noise terms as 
well, we modify our unitary control operator to 



(20) 



where Uz is another non-zero integer satisfying the condi- 
tion Tlx ^ ?T-z • The unitary control operator now consists 



(k) 

of two rotation operators. The cri part of the unitary 
operator rotates the cri'^'' and Uy'^'' noise terms and aver- 

(k) 

ages them out to zero, while the tri part of the unitary 
control operator takes care of the remaining ci'^'' noise 
terms (it nevertheless also rotates the cry*'' operators). 
The condition Ux 7^ nz is important because otherwise, 
the effect of the second rotation cancels some effect of 
the first rotation such that, for instance, the ay''^ noise 
terms do not average out to zero. All these claims can 
be examined by explicitly verifying if Eq. ([5]) holds. For 
instance, we observe that 



(21) 



/o 
Also, 



tc 



C/t(i)crW[/^(i)dt 



2irm^£7<'°'t/tc„47rmjccr^'°'t/tc „-27rin^(T<'°'t/tc„(fc) jj. 

O fc- fc' tit/. 







(22) 



Using the commutator 

[g2«..„.af)t^g-»..n.4'"tj ^ _2isin(2n,wi)sin(n,wt)4'^), 
where a; = ^fJQ can simplify, 

Ul{t)(T^^^Uc{t)dt 



(23) 



which further simplifies to 



2i / [cos{nzUJt) + ism{nzLi!t)(T 



(fc)l 



Now, 



X [sm{2nxUjt)sm{nzU}t)]al'''>ai''Ut. (24) 



ain^ (nzUit) sm(2nxUJt)dt = 0, (25) 



no matter what the values of Ux and Hz are, but in order 
to have 



cosluzUjt) sm{2nxi-ut) siii{nzUjt)dt = 0, (26) 



we require that Hx ^ tiz- Therefore, if Ux tIz^ the 
noise terms are eliminated. Similarly, one can check that 
the cTy'^' noise terms are also eliminated. 
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The necessary control field to implement the unitary 
control operator Uc{t) can be found from the Schrodinger 
equation = Hcit)Uc{t). We then find 

2 



i=l 

con. 



(27) 

One obvious aspect of this control field is that both qubits 
are addressed in exactly the same way. The field config- 
uration is also quite simple: it consists of a local static 
field and a local rotating field. 

However, the control Hamiltonian found above is not 
universal. In particular, it cannot eliminate all possible 
forms of system-environment coupling shown in Eq. ([2]). 
For instance, consider the noise term proportional to 

cri^Vi^^ We find that 



g47rm, fT_li) t/t, g4vr™. o-i" t/to ^(2) 

3s{2nzUJt) + i sin(2n3wt)cr^2) a^^^a^J^Ut ^ 0. 

(28) 



This is obvious because 



cos 



■{2n,ujt)ai^')cri^^dt ^ 0. 







Therefore, the ui^Vi^' noise term does not average out 
to zero if ai^^ and are rotated at the same frequency. 



B. Universal protection of two-qubit states 

We have just shown that it is not possible to eliminate 
all the noise terms by applying the same field to both 
qubits. So the important question is the following: is it 
possible to find a field configuration in which, by applying 
different fields to the two qubits, all the noise terms as 
shown in Eq. ^ can be eliminated? At the same time, we 
would also like to retain the relative simplicity of the field 
configuration used for previous local noise considerations. 
This motivates us to investigate if 



us)^u^'Ht)uPit), 



(29) 



where 



serves to eliminate all the noise terms. Note that, since 
we allow the possibility of different fields being applied 



to the two qubits, this Uc{t) differs from the previous 



Uc{t) in that, this time, the integers in Uc^^ and [/, 



r(2) 



need not be the same (previously we had niP = n)^' and 
n'z'^ = ni^^). The postulated Uc{t) is obviously periodic 
in time with period ic- Furthermore, as shown below, we 
find that all the noise terms can indeed be eliminated, 



,(2) 



provided that the integers n>x \ n^z \ n^x\ 
some criteria. For simplicity, we consider n'"x \ n,i^\ nx\ 
and riz to be positive integers. Also, since we expect 
that the integers are different, for our own convenience 
in narrowing down the criteria fulfilled by the integers, 
we impose the condition 



and n^z^ fulfill 



,(1) 



< nl^^ < n 



(2) 



< n 



(2) 



(31) 



Let us now find the criteria that the integers n"x \ n^z\ 

(2) (2) 

Tlx , and Tiz need to fulfill. In order to do this rigorously, 
we need to check that each noise term averages out to 
zero under the action of the applied fields. Since we have 
ordered the integers as in Eq. (PT|) . we already have that, 



(1) 



(2) 



(32) 



(33) 



Using the derivations in Sec. II-A, it is easy to see that 
all local noise terms are indeed eliminated. 

We next examine the fate of nonlocal noise terms. For 
instance, let us consider the noise term proportional to 
cf^x'' cf'x^ ■ This time we have, 



Ul{t)a^}^a'i^U,{t)dt 

g47rmii'cr(i)t/t,^(l)g47rin(,2't7(='t/t,^(2)^^ 

)s(24i)wt) + is\n{2n'^^'^ujt)a^p 
cos(2n(2'a;t) + i sin{2n'f'^ ujt)a'^z^ 



which is zero, provided that, 



(34) 



(35) 



This is obvious because under the condition n'z^ ^ n'z\ 



i{2n^z^(^jt) cos{2n^z^ uj't) dt 



and similar terms are all zero. 

Therefore, one observes that if cri^'' and are rotated 
at different frequencies, then the noise term proportional 
to cri^'cri^'' is eliminated. This condition gives support to 
our intuition that the fields applied to each qubit should 
be different. 
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We next outline the calculation for the noise term pro- 
portional to . In this case, the calculation is con- 



unitary control operator in Eq. ([29]) with the following 
conditions, 



sidcrably more involved. To calculate the required in- 
tegral involving this noise term, we first calculate (sup- 


< < 




pressing the k index). 






4^) 


+ 










4^) 


+ 










4^) 


+ 










4^) 


+ 




+ sin(2n:ca;t) cos(2nza;i)crj,. (36) 






— n 


(2) 

Z 




The integral that wc wish to set to zero then becomes 






+ n 


(2) 


^0, 



[/t(t)4i)42)t/^(i)dt 

cos(24^'w^)4^' + sin(24^'wt) sin(24^^wt)4^^ 
+ sin(24^^wt) cos{2n'^pujt)a''y^ 

i{2n'^^ujt)af^ + sin(24^^t^i) SYa{2n'i^ ujt)af'i 



+ sin(24^^w0 cos(2n(2)cjt)CT(2) 



(37) 



By multiplying the terms in the square brackets above, 
we get different terms. Each of these terms must indi- 
vidually integrate to zero, because the tensor products of 
two Pauli matrices are linearly independent in the oper- 
ator space. So, for example, wc require that 

cos(24^^wO sin(24^^t^i) sin(24^'wt)4^Vi^^dt = 0, 

(38) 

which is true provided that 

4^)-t-4^)-4^)^o. (39) 

One might think that we would also need three other 
conditions, one of which is given by 



,(1) 



42) + n^^^ ^ 0. 



(2) 



(40) 



Fortunately, due to the ordering to the integers in 
Eq. ([5T|) . this condition and the other two are redundant. 
Therefore, we can ignore these redundant conditions. 

Similarly, analyzing each of the others terms in Eq. ([37]) 
one by one, and keeping the ordering of the integers in 
mind, we arrive at the following list of criteria: 



4^)-!- 4^) 



4^) ^4^) + 4^), 
4^^^4^^ + "W, 
4^^ - 4^^ ^ 0, 

,(2) ^„(2 



,(1) 



,(1) 



- nl 



^0. 



(41) 



The other seven types of system-environment coupling 
shown in Eq. ([5]) can be treated in a similar fashion and 
will not be repeated here. Carefully going through all of 
them, we come to the conclusion that by applying the 



(42) 

our two-qubit system can be (approximately) decoupled 
from the most general environment. From the above con- 
ditions, it is seen that not only must the frequencies in 
Uc{t) be all different, but also that neither of the two 
larger frequencies should be the sum of two smaller fre- 
quencies. Furthermore, the difference of the two larger 
frequencies should not be equal to the sum or the differ- 
ence of the two smaller frequencies. It is not hard to find 
integers that fulfill all the criteria we have found. One 
possible choice is n^x^ = 1, ni^^ = 2, nl^^ = 4, ni'^' = 8. 

Finally, the control Hamiltonian, Hc{t), which is 
needed to generate the unitary operator Uc{t), is found 
to be [from the time derivative of Uc{t)] 



-|- uin 



[cos(2lj4^04*^ - sin(2a;4^04'^] } • 

(43) 



Each of the two qubits is now subject to a different lo- 
cal control field consisting of a static field and a rotating 
field - we must address each qubit individually. With 
these control fields, the two-qubit system is dynamically 
decoupled from the environment, for all possible types of 
system-environment coupling. Note also from Eq. (|43p 
that the field amplitude should also go up if the frequen- 
cies of the driving field are increased to compete with the 
cutoff frequency of the environment. 



III. PROTECTION OF TWO-QUBIT GATES 

Once we have the control operator Uc{t) that is able 
to protect two qubits against decoherence in a univer- 
sal manner, the next natural question is how to turn on 
coherent evolution in two-qubit systems such that two- 
qubit gates can be also protected. This is important 
because, in reality, there is no instantaneous quantum 
gate. As shown below, we can extend our previous con- 
siderations to protect a two-qubit state and implement a 
desired gate at the same time. Some early studies con- 
sidered the protection of a two-qubit gate against ran- 
dom dcphasing 52 1 and against bit-fiip errors [s^, but 
these early decoherence suppression approaches are not 
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applicable to an arbitrary environment. Our procedure 
is analogous to Ref . H^, |43| , but for most general system- 
environment coupling in two-qubit systems. The exten- 
sion here is worthwhile because in actual realizations of 
two-qubit gates, it is unavoidable that the two-qubit in- 
teraction Hamiltonian will suffer from fluctuations, on 
top of local noise terms seen by each individual qubit. 



A. Two-qubit gate under pure dephasing 

To illustrate the method, we start off with the simple 
case of pure dephasing. The interaction between the two 
qubits and their environment is given by. 



H. 



dephasing 



B^)a^+Bf)af)+B^)a^af^. (44) 



As can be easily verified, in this simpler control 

operator 

Uc{t) ^ cx-p{~2T:ia'^^nit/tc) cxp{~2T:ial^^n2t/tc), (45) 

with rii ^ suffices to protect two-qubit states. 

Consider now a two-qubit gate that converts a separa- 
ble state into a Bell state, i.e.. 



where 



1 



V2 



(It). + 1;> J i;>. 



f 



V2 



(in). + i;t) J, 



In order to find this new Hamiltonian, we begin by 
writing the system Hamiltonian as 



(49) 



where Ho(t) implements the gate. In order to find the 
unitary control operator that both implements the gate 
and protects against decoherence, the basic idea is to 
once again transform to the frame given by Hc{t). Now, 
in this frame, the effect of the environment has already 
been largely removed - it is almost as if the environment 
were not there. Therefore, we implement the gate in this 
picture. After doing so, we simply transform back to our 
original reference frame to find the total unitary control 
operator. 

Let us now carry out these ideas in detail in order to 
find the required Hs{t). We first, once again, write the 
total Hamiltonian as 



(50) 



where H'{t) = HQ{t) + H-q + i?sB- We now transform 
to the frame of the control fields, as we did before. In 
this frame, Ha{t) becomes Ho = Ul{t)Ho{t)Uc{t). Cor- 
responding to this Hamiltonian, there is a unitary time- 
evolution operator. 



Ua{t)=T exp 



Ho{s)d.8 



(51) 



lt). = ^(|0) + |l)), 

li>. = ^(|o)-|i)), 



(46) 
(47) 



with |0) and |1) being eigenstates of the operator. 
We consider the initial state to be \ipo) in order to bring 
out the effect of the dephasing noise clearly. It should 
be noted that the above gate is analogous to the usual 
controlled-NOT (CNOT) gate, since the CNOT gate per- 
forms the operation -ij(|0)-f |1))|1) — > ^(|01) |10)). 
Therefore, we refer to the gate implementing the above 
operation as the CNOT gate. We work with the CNOT 
gate because it generates entanglement - the usual CNOT 
gate acting on -^(|t).+li).) \i)x yields a separable state. 

The most straightforward way to implement the 
CNOT gate (up to an irrelevant global phase) is to use 
the two-qubit Hamiltonian, 



2t 2 ' ^ 



(48) 



where r is the time over which the gate is implemented. 
Note that no decoherence control fields are being applied 
at this stage. Therefore, during the gate operation time, 
the two-qubit state is vulnerable to decoherence due to 
the environment. Our task is to modify the Hamilto- 
nian given by Eq. (|48p such that the new Hamiltonian 



not only implements the CNOT gate, but also prevents 
decoherence. 



It is this unitary operator that we use to implement the 
gate. Therefore, 



Uo(t) = exp 



, (52) 



where again r is the time over which the gate is imple- 
mented up to a global phase, that is, C/gato ~ Uo{t ~ t) . 
We set T = Ntc {N is a positive integer). Comparing 
Eqs. dST]) and it is clear that 



TT 1 



(1) 



r(2) 



(53) 



will do the right job (this choice for Hq is simple because 
it is time- independent). 

As shown in our previous general consideration of DD 
in Sec. II, the total system-environment time-evolution 
operator in the frame of the control fields is already ap- 
proximately decoupled into a product of system and en- 
vironment parts [see Eq. (jl5|) ]. In particular, applying 
the Magnus expansion to the following total evolution 
operator 



Texp 



-I I H'{s)ds 



we have that for sufficiently small tc, 



-iHsT 



(54) 



(55) 
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Wc finally transform back to the original frame. In this 
frame, the unitary time evolution operator is given by 



C/tot(T) « C/c(r)[/gatce- 



But Uc{t) is just identity, leading to 



gatc^ 



(56) 



(57) 



Clearly then, the desired gate operation is performed on 
the two-qubit system. 

For arbitrary time t, the system is also approximately 
decoupled from the environment (because in the limit 
tc « 0, t is always close to an integer multiple of tc). 
Then, in the lab frame the overall unitary evolution op- 
erator for the two-qubit system at arbitrary time t is 
given by 



Usit) = Uc(t)Uo{t) 



= exp{-2TTia^J'nit/tc) cxp{~2TTi(T^J>n2t/tc)Uo{t). 

(58) 

Further using the time-dependent Schrodinger equation, 
the Hamiltonian that generates the overall evolution op- 
erator Us{t) can be obtained as follows. 



Hs{t) ^ ujmai^^ 



TT 1 - 

cr^^^ cos(2wn2t) - cr^^^ sin(2a;n2i) 
cr(,^Vf ) cos(2wn2<) + cri^V^^' sin(2wn2t) 



(59) 



with ni 7^ 712. By our construction above, such a field 
configuration implements the gate and protects against 
two-qubit pure-dephasing at the same time. Note that 
here some nonlocal field components are needed. This 
is expected. After all, the original CNOT gate Hamilto- 
nian, Eq. also needs a qubit-qubit interaction term. 
The message is that an oscillating qubit-qubit interaction 
can be highly useful in implementing robust two-qubit 
gates in a noisy environment. This need for oscillating 
qubit-qubit interaction here should not be regarded as 
a great disadvantage of our universal continuous DD. In 
fact, even in pulsed DD schemes for entanglement pro- 
tection pulsed qubit-qubit interaction is necessary 
to reduce the number of UDD layers. The requirement 
for time-dependent qubit-qubit interaction terms is also 
consistent with previous case studies under the general 
DCG framework [ill, El. 



B. CZ gate protected against most general 
environment 

We are now ready to carry out similar calculations to 
construct a field configuration that protects a two-qubit 
gate against all possible forms of system-environment 
coupling. As an example, we consider the implemen- 
tation of a controlled phase (CZ) gate [SJ]. This case is 



representative because, if we can reliably implement the 
CZ gate in the presence of arbitrary decoherence sources, 
then, together with single-qubit gates, we can perform 
universal gate operations in the presence of an unknown 
environment. 

The CZ gate, in a matrix form in the standard repre- 
sentation, can be written as, 



Ucz = 



/ 1 
10 
10 

Vo -1 



(60) 



In order to achieve this unitary operation (up to a global 
phase) on a two-qubit system, a desired bare Hamiltonian 
would be. 



TT 1 

2^2 



(61) 



where t is the gate operation time. But the resulting two- 
qubit state is not protected against the environment. As 
such, we seek instead a time-dependent system Hamilto- 
nian Hs(t). 

Our previous treatment for the CNOT gate in a pure 
dephasing model can be extended easily. The physi- 
cal picture underlying the technique remains exactly the 
same. That is, we implement the desired gate in the ro- 
tating frame and then transform it back to the lab frame. 
Following our previous notation, in the rotating frame we 
hope to have. 



Uoit) = exp 



. TT i 

"272^' 



r(l) 



.(2) 



.(1)^(2) 



(62) 



On the other hand, Uc{t) is given by Eq. ([29| . By com- 
bining these two unitary operators as before, the sought 
Hamiltonian is determined by simply using the time- 
dependent Schrodinger equation on the overall unitary 
evolution operator. It is found to be, 



H^{t) = {w^'^f^x^ [cos(2u;n«t)a« 

TT 1 



i=l 

sin( 



in(2c.r.Wt)a«] } + ^- [ai'^ cos(2^4i)t) 



,(1) 



(2), 



,(2)y 



- cri^V(2) cos{2ujn2-h) cosi2ujn'^^h) 
+ cri^V(2) cos(2w4^)t) sin(2a;4^)t) 
+ crt^Vi^) sm(2LJni^h) cos{2Ljn^^h) 

- smi2ujni^h)sm{2un<;^h) 



(63) 



We stress that here we did not make any assumption 
about the system-environment coupling. A CZ gate can 
hence be implemented and protected against any type 
of decoherence, so long as the driving frequencies are 
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sufficiently large (also sufficiently strong) relative to the 
cutoff frequency of the environment. Comparing the sys- 
tem Hamiltonian here with that in the previous pure- 
dephasing case, more oscillating qubit-qubit interaction 
terms are required for decoherence suppression. Another 
interesting observation is that here, the oscillating qubit- 
qubit interaction terms carry the sum and the difference 
frequencies 2{n^x'' + n>x^)uj and 2(ni"'^' — ri"x^)uj. This 
feature can be regarded as a result of the dual role of 
the control fields (implementing and protecting a gate). 
It is also consistent with the fact that the two qubits 
should be rotated at different frequencies. We do not 
suggest that the required control fields in Eq. (|63| are 
easy to realize experimentally. But at least, such an ex- 
plicit solution as an example of universal DD is indica- 
tive of what could be crucial in protecting two-qubit gates 
without making assumptions of system-environment cou- 
pling. Our two-qubit gate construction also constitutes 
an explicit and simple implementation of DCG [4l|, |42| 
to fight against arbitrary (environment-induced) single- 
qubit and two-qubit errors, using a static field plus sev- 
eral continuous-wave driving fields of different frequen- 
cies. 



IV. THE MASTER EQUATION 

In this section, we bricfiy introduce a master equation 
for our use in numerical simulations. Here we only sum- 
marize the basic formalism. For a conceptually simple 
derivation of the master equation based solely on time- 
dependent perturbation theory and not on any special 
techniques such as projection-operator methods, please 
see Appendix B. 

In the master equation approach [ssj , we find the dy- 
namics of a system that is coupled to the environment by 
first considering the total system consisting of the system 
and the environment as closed. We then trace over the 
environment to obtain a differential equation, known as 
the master equation, for the reduced density matrix of 
the system. 

We start by writing down the Hamiltonian of the total 
system as 



Htot = Hsit) + Hb + HsBit), 



(64) 



where Hs{t) describes the system Hamiltonian, Hb the 
environment Hamiltonian, and Hqb is the interaction be- 
tween the system and the environment. In general, we 
can write 



i/sB(t) = $^i^,(0®S.(i)- 



(65) 



Here, the Fj{t) are operators in the system Hilbert space, 
and the Bj(t) are operators in the environment Hilbert 
space. 

We assume that the interaction between the system 
and the environment is weak, and that the total initial 



state of the system and environment is a product state. 
We can then derive the master equation describing the 
time evolution of the reduced density matrix p of the 
system as. 



dpjt) 
dt 



iW),H{t)] + 



f dsm{t,s)pit),F,it)]Ci + h.c.}, (66) 



where, 



F,(t,s)^Usit,s)F,is)Ul{t,s), 



Cl^{B,{t)B,{s)), 



B,{t) = ul{t,to)B{t)UB{t,to), 



(67) 



(68) 



(69) 



and UBit, to) and Us{t, to) are the unitary time-evolution 
operators corresponding to Hb and Hs{t) respectively. 
Such a master equation has been used previously in Rcf. 

M- 

We consider the environment as a collection of an in- 
finite number of harmonic oscillators, so that 



(70) 



Here index k denotes different modes of the oscillators in 
one bath, and index j denotes different thermal baths. 
Furthermore, we take the Bj operators for the jth bath 
as 



(71) 



where the gj^k are coupling strength parameters. All the 
baths are assumed to be in a thermal equilibrium state 
with the same temperature T. We then note that the 
bath correlation function, given by Cf^, can be written 
as a function of the time difference t — s. In order to pro- 
ceed with the calculation the bath correlation functions, 
the discrete modes of the environment are replaced by 
a smooth continuum of modes specified by the so-called 
spectral density J(f2). For our numerical simulations, we 
consider an Ohmic spectral density with an exponential 
cutoff, that is. 



j{n) = one 



(72) 



where G is the coupling constant, and Wc is the cutoff 
frequency of a bath. For simplicity, we assume that all 
baths have the same spectral density, with the same G 
and ujr- 
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NUMERICAL RESULTS 



In this section wc present numerical results illustrating 
universal protection of two-qubit states and two-qubit 
gates. Two extreme cases are considered. In the first 
case, two qubits are coupled to fifteen different baths, 
each of which induces one type of system-environment 
coupling. The operators Fj in the master equation (|66|) 



arc then given by. 



0) 

o 



3 
o 

o 
O 




(73) 



where ctq = /, and the notation for Pauli matrices is the 
same as before. In the second case, all possible types of 
decoherence are modeled by a common bath. That is, in 
the master equation [see Eq. (|66p ]. there is only one Fj, 
which is written as. 



> k=l 



.(1) 



3 

k=l 



A2) 



3 

E 

k=i,i=i 



(74) 



In addition, pure-dephasing cases arc also considered, 
with the different and common bath cases defined in 
an analogous way. The only difference is that for pure- 
dephasing cases, there are no terms containing or ay 
in the Fj operators. 



Wc work in dimensionlcss units with h = 1 and 
ksT = 2. In these units, the parameters we use are 
LUc = 27r, and tc = 0.5 (so oj = 47r), unless stated 
otherwise. As the measure of bipartite entanglement, 
we use the concurrence [57| . Given a two-qubit den- 
sity matrix p, the concurrence, C, is defined as C = 
max{Ai — A2 — A3 — A4,0}, where the Xi are the square 
roots of the eigenvalues (in descending order) of the ma- 



trix p{all^ (g) <^'y')p*{oy''' ® <^'v') (the asterisk denotes 
complex conjugation). 

We first present results of two-qubit state protection 
using our universal continuous DD fields. For conve- 
nience, the self-Hamiltonian of the two-qubit system is 
set to zero. 



(2)^ 



.(1) 



2)^ 



FIG. 1. Entanglement vs time (in dimensionless units) with- 
out decoherence control fields (dotted line) and with applied 
control fields (dashed and sohd fines). The environment is 
modeled by 15 different baths, i.e., in our master equation, 
the system coupling operators are given by Eq. (|73p . The 



dashed line is for rix 



(1) 



1, ri; 



(1) 



2,n: 



(2) 



(2) 



8. 



The solid line is for stronger and higher frequency fields with 

16. For this numerical 



ni*' = 2,ni^' = 4,71, 
simulation, we use G 



(2) 



,(2) 



0.05. 



Figures [T] and [2] illustrate the performance of the de- 
coherence control fields in protecting two-qubit entangle- 
ment against an environment that generates all types of 
decoherence. Without these fields, we see (dotted curve) 
that, in both the common-bath and different-bath cases, 
the entanglement rapidly decays due to the interaction 
with the environment. The situation changes dramati- 
cally after switching on the continuous control fields. The 
dashed curves demonstrate the suppression of entangle- 
ment decay due to the control fields. Furthermore, by 
applying fields of greater strength and higher frequency 
(the solid curves) - thus effectively reducing tc - even 
better protection of entanglement is achieved. 



CD 

o 

<D 
i_ 
i_ 

o 
c 
o 
O 




FIG. 2. Entanglement versus time (in dimensionless units) 
without decoherence control fields (dotted line) and with ap- 
plied control fields (dashed and solid lines). Here the envi- 
ronment is modeled by one common bath, with the system 
coupling operator given by Eq. (|74[) . The parameters used 
are the same as in Fig. [1] 
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We now study the effectiveness of the control Hamilto- 
nian found in Eq. (|59p in implementing the CNOT gate in 
the presence of pure dephasing. First of all, the Hamilto- 
nian that only implements the gate without decoherence 
control is given by Eq. (j48p , which we rewrite here for 
convenience, 

This Hamiltonian should be contrasted with the Hamil- 
tonian in Eq. ([5^ that both implements the gate and pro- 
tects against decoherence. Using numerical simulations, 
we can directly compare the performance of these two 
Hamiltonians for entanglement generation in the pres- 
ence of pure dephasing. 
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FIG. 3. For the CNOT gate, behavior of entanglement versus 
time (in dimensionless units) using a bare Hamiltonian that 
only implements the gate (dashed line) and using control fields 
that both implement the gate and protect against all types of 
pure-dephasing (solid line) for the case of different baths. In 
the dimensionless units defined before, the parameters used 
are m = 2, n2 = 1, r = 0.5, G = 0.03. 
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FIG. 4. Performance of the CNOT gate in generating entan- 
glement using a Hamiltonian that only implements the gate 
(dashed line) and using control fields that both implement 
the gate and protect against all types of pure-dephasing at 
the same time (solid line) for the common bath case. The 
parameters used are the same as for Fig. [3l 



Such a comparison is done in Figs. [3] and |4l The 
dashed curves depict the entanglement generation by 
the bare Hamiltonian Eq. (|48|) and the solid lines are 
for the performance by the control Hamiltonian given 
in Eq. ((59)) . As expected from an entangling gate, both 
Hamiltonians generate entanglement with similar per- 
formance in the beginning. However, after some time 
(for common-bath and different-bath cases), the effect 
of the environment is accumulated and eventually the 
bare Hamiltonian loses its battle against the environ- 
ment, whereas for our constructed control Hamiltonian, 
the entanglement generation stays close to its expected 
value. By the time the gate operation is completed, much 
better performance is achieved due to the application of 
continuous DD fields. 

We now perform a similar task for the CZ gate in the 
presence of all types of decoherence. The CZ gate is 
used to take a separable state to a fully entangled state. 
We compare the gate performance afforded by the con- 
trol Hamiltonian given by Eq. (|63|) with that of a bare 
Hamiltonian that solely implements the CZ gate. The 
bare Hamiltonian, as already given by Eq. (|6ip , is 

2t 2 V - z z ) 

Up to a global phase factor, this Hamiltonian implements 
the CZ gate in time r. We stress that numerical simu- 
lations here are no longer restricted to pure dephasing. 
Instead we are considering the most general case, allow- 
ing errors such as uncorrelated bit flipping and dephasing 
as well as 'noisy' interaction between the two qubits. 
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FIG. 5. Evolution of entanglement using a bare Hamiltonian 
that only implements the gate (dashed line) and using con- 
trol fields that both implement the gate and protect against 
decoherence at the same time (solid line) for the case of 15 
different baths [Eq. (|73|) ]. The parameters (in dimensionless 
units defined before) are r = 0.5, G — 0.02, ujc = tt, ni^' = 
l,7ii^' = 2, Tii^' = 4, nl^' = 8. The CZ gate here converts a 
separable state to an entangled state. 

In Figs. [S] and [Bl the dashed curves depict the per- 
formance of the bare Hamiltonian in Eq. (j6ip , whereas 
the solid curves show the performance of the control 
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Hamiltonian we found in Eq. (|63|) . Once again, the per- 
formance benefit is obvious. In both the different-bath 
and common-bath cases, with the continuous DD fields 
implemented, we are able to achieve almost perfectly en- 
tangled states even in the presence of all possible types 
of decoherence. By contrast, the desired coherent evolu- 
tion takes place with a clearly poor fidelity if only a bare 
Hamiltonian is used. 

1 I 1 1 1 1 — ■ 
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FIG. 6. Considering the CZ gate, evolution of entanglement 
using a Hamiltonian that only implements the gate (dashed 
line) and using fields that both implement the gate and pro- 
tect against decoherence at the same time (solid line) for the 
case of common bath [Eq. (|74|l ]. The parameters used are the 
same as in Fig. (5] 



VI. CONCLUSION 

In this paper, we first asked and answered the following 
question: is it possible to use continuous fields to achieve 
(albeit low-order) protection of two-qubit states as a uni- 
versal dynamical decoupling approach? By extending the 
methodology in Ref. [49[, we have found a rather simple 
field configuration to achieve this task. This associated 
decoherence control is completely general in the sense 
that it is able to protect the state against all types of de- 
coherence, so long as the frequency of the control fields is 
sufficiently large (with sufficient field strength) as com- 
pared with the environment cutoff frequency. From a 
practical point of view, the very existence of a universal 
scheme is important if we do not have enough informa- 
tion about the environment. The found continuous DD 
is also relatively simple - only local continuous and pe- 
riodic fields are required. Our results are thus at least 
complementary to recent studies of universal pulsed DD 
for entanglement protection. In particular, under the 
circumstances where multi-pulse DD is difficult to im- 
plement (e.g., due to the requirement of very small pulse 
intervals), then our universal continuous DD scheme pro- 
vides one alternative. Furthermore, one can imagine us- 
ing a combination of pulsed DD and continuous DD to 
reliably store quantum information. 

We have also constructed continuous control fields to 



implement two-qubit gates in the presence of most gen- 
eral decoherence. This is important for three reasons. 
First, it always takes a finite amount of time for a quan- 
tum gate to be implemented and as such a two-qubit 
gate must be protected against decoherence during the 
gate operation time. Second, during the implementation 
of a two-qubit gate, coherent evolution of the system it- 
self complicates the issue of decoherence control due to 
the transformation between different types of quantum 
coherence properties. Third, the implementation of two- 
qubit gates itself will unavoidably bring about noise in 
qubit-qubit interactions. As seen from two case studies 
of universal two-qubit gate protection against both lo- 
cal and nonlocal noise, the required control Hamiltonian 
is not too complicated, with the most involving compo- 
nent being oscillating qubit-qubit interaction terms. Our 
treatment is general in the sense that we have not con- 
sidered any particular physical implementation of a two- 
qubit gate. It would be interesting to apply our find- 
ings here to a particular physical realization of two-qubit 
gates. One excellent example would be in the recent im- 
plementation of superconducting two-qubit gates using 
simple microwave fields [ssj . In a second example from 
the trapped- ion context, continuous microwave driving 
is already theoretically considered to protect two-qubit 
gates against noise due to magnetic field fluctuations and 
the thermal motion of the ions [H^- Of course, in such 
physical realizations, it may be the case that only a few 
noise sources contribute appreciably to decoherence and 
therefore the required continuous DD fields may be sim- 
plified. Our results here also lay a useful starting point 
for future optimization studies [13, HH, by, for exam- 
ple, first extracting some information about an environ- 
ment. Finally, as pointed out earlier, this work offers an 
explicit and simple route to construct dynamically cor- 
rected two-qubit gates [4l], \^ to fight against arbitrary 
system-environment coupling. 
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Appendix A: The Magnus Expansion 

The Magnus expansion (sH says that the unitary evo- 
lution operator, U{t), corresponding to a time-dependent 
Hamiltonian, H(t), can be written as 

oo 

[/(t) = exp^^,(i) (Al) 
1=1 
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where 



Ai = -i [ dtiH{ti), (A2) 

"'0 

A2 = -l f dt, f ' dt2[H{t,),H{t2)], (A3) 
^ Jo Jo 

with higher order terms given by higher order eommuta- 
tor expressions. For further details, we refer the reader 
toRef. [62|. 



Appendix B: Derivation of the master equation 

Here we present a conceptually simple derivation of the 
non-Mar kovian master equation that we have used. Such 
a master equation has been previously used in, for exam- 
ple, Ref. [53. However, unlike Ref. [5g, we do not use 
any advanced techniques such as the projection-operator 
method to derive the master equation. Rather, we only 
use a simple perturbation theory. Note that it is essential 
that we do not make the Markov approximation since we 
are interested in decoherence control. 

For simplicity, here we derive the master equation for 
a system interacting with an environment consisting of a 
single bath. The more general case of the environment 
consisting of multiple baths can be dealt with via a simple 
extension. As usual, we start off by writing down the 
Hamiltonian of the total system. 



Htotit) = Hit) + V{t). 



(Bl) 



H{t) ~ Hc,{t) + is the Hamiltonian describing the free 
system and the bath, and, for notational convenience, 
here we use V{t) instead of Hsb to represent coupling 
between them. Note that these are Schrodinger picture 
operators. Any time dependence in the system Hamilto- 
nian or the Hamiltonian describing the coupling is an ex- 
plicit time dependence. This time dependence can arise, 
for instance, due to the application of external fields. 

Consider system-environment operators of the form 
Y ^ Ib, where Y is an operator acting on the Hilbert 
space of the system and Ib denotes identity in the Hilbert 
space of the environment. If the state of the total sys- 
tem is described by the density matrix ptotit), then the 
expectation value of the operator Y is given by 



(y) =Trs(rp), 



(B2) 



where p = Treptot- Ttb denotes taking trace over the 
bath degrees of freedom, while Trg means tracing out 
the system degrees of freedom. Now, p is our primary 
object of interest. We are trying to derive an equation of 
motion for p. Note also that 



(Y) -Trs3[(r®/B)/9tot]. 



(B3) 



Now, for computational purposes, we express the density 
matrix in some basis, namely, Pmn{t) = (™| p(0 I"-)- 



particular, if we choose Y = \n) {m\ = Ynnn we get 

(y) = Trs[r„,„p(t)] 

= (m| p{t) \n) = pmn{t) 

= Trs,B[(F„„(8)/B)ptot(t)] 

= Trs,B[f/^(t,io)(Km <»/B)C/(t,to)ptot(io)], 

(B4) 

where U{t^to) is the unitary evolution operator describ- 
ing the unitary evolution of the total system. The 
cyclic invariance property of the trace has also been 



used. Defining F„. 



Ib 



Xr, 



we observe that 



U^{t,to)XnmU{t,to) is just & Heiscuberg picture oper- 
ator. We refer to this operator as Xnm{t), with the un- 
derstanding that any X operator with a time dependence 
is in the Heisenberg picture. Therefore, we can write. 



Pmn{t) = Trs3[ptot(io)^' 
It follows that 

dpvin{t) „ r ,. .dX 

T7 



dt 



Trs.i 



Ptot {to 



dt 



(B5) 



(B6) 



We now derive the Heisenberg equation of motion for 
a general system-environment operator, and substitute 
that in the above expression. After doing so, we take 
the trace over the system and the environment to obtain 
our master equation. We proceed with the derivation of 
this Heisenberg equation of motion using perturbation 
theory. 

As in standard time-dependent perturbation theory, 
we set U{t,to) = Uo{t,to)Ui{t,to) where Uo{t,to) = 
Us{t,to)UB{t,to) describes the free evolution of the sys- 
tem and environment, and the Ui{t,to) comes in due to 
the coupling with the environment. It follows then, that 
to first order in the system-environment coupling. 



(B7) 



U{t,to)^Uoit,to) [l-i dsVis) 



to 



where Vis) = U^ois,to)Vis)Uois,to). 

We know that operators in the Heisenberg picture and 
operators in the Schrodinger picture are related by. 



OHit) = uHt,to)Oit)U{t,to), 



(B8) 



where 0(t) is a general Schrodinger picture operator, and 
Onit) is its Heisenberg picture counterpart. 

Considering Oh {t) to be a general system-environment 
operator that has no explicit time dependence, we have 
the Heisenberg equation of motion. 



dOHjt) 
dt 



l[HH{t),OHit)]+l[VH{t),OHit)]. (B9) 



Using our above perturbative expression for the uni- 
tary time evolution operator, and observing that 
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[Vff (t), Onit)] is a Heisenberg picture operator, we have, 

dOHit) 



dt 



i[HH{t),OH{t)]+i[V{t),0{t)] 
+ f ds[[V{t),6{t)lV{s)]. 

J to 



(BIO) 



We now set Onit) = Xnm{t), and substitute Eq. (|B10[) 
in Eq. (|B6p . Each term of Eq. (|B10p gives a term in the 
master equation, so we work them out one by one. We 
assume that the initial state is a product state p(to)® Pb- 
The first term of Eq. (|BT0)) . i.e., i[HH{t),Xnrn{t)], 
when substituted in Eq. (|B6|) . leads to 



Trs,B{(p(io) ® PB)i[HH{t),Xnrn{t)]] 
— ^ ^ ^ ^ra7ira' 7i' Pra' n' -i (-^-^-^) 

where /^.^fuiYn' n' = ^mm'-^n'n{i^ ^rui' -^7nrn' The SeC- 

ond term of (jBlOl) . i.e., i[V"(t), X„m(t)], leads to 

Trs,B{*(p(to) ® /OB)f/o^(t, toWit), Xnm]Uo{t, to)}. 

The system-environment coupling is of the form, 

V{t) = F{t) (g) B{t), (B12) 

where F{t) is an operator acting in the Hilbert space of 
the system and B{t) is an operator acting in the Hilbert 
space of the environment. With such a coupling, we can 
then work out the trace over the environment. This is 
given by 

Ti-B[pBUl{t,to)B{t)UB{t,to)] = {B{t)). (B13) 

We assume that this is zero. This is commonly referred 
to as 'centering' of the bath. 

The last term of (|B10p . when substituted in (jB6[) . gives 
us four terms. Here, we only show the working for one 
of them. The rest can be worked out in almost the same 
way. We now proceed to simplify 



Trs3[(p(to)®PB) / dsV{t)X„^{t)V{s)]. (B14) 

Jto 

The trace over the environment gives 

Tre [pBUlit, to)B{t)UB{t, to)Ul{s, to)B{s)UB{s, to)] 
= {B{t)B{s)) =Ct,. (B15) 

The trace over the system gives 

Trs[pit)F{t)YnrnUs{t, s)F{s)Ul{t, s)]. 

We can now use the completeness relations of the basis 
states to write 



-F(,^) ^ ^ -^n' 7n' it)Yrfii j^i' ^ 



(B16) 



F{s)=Y,F^0{s)Y^0, 

af) 



(B17) 



(B18) 



Ul{t, s) = Y. uI'-'^'y,... = C/r'"V,v'. (B19) 



Substituting these relations back, and assuming that our 
basis states are orthonormal, we can simplify the trace 
over the system to 

^ ^ K'„(i)^^™>' (.s)C/F"' (i, s)vY''' (t, S)p,'n' (t), 

n'm' 

where p{t) has been replaced by p{t). This is justified 
since the correction gives us terms of higher order in the 
coupling strength in the master equation. For conve- 
nience, we define, 

gX,it, s) EE s)uf^'it, s). (B20) 



Putting all the above expressions together, Eq. (|B14p 
simplifies to 

7l'm' l-LU 

(B21) 

After working out all other terms, we finally arrive at the 
master equation 



" ^ ^ ^ ^rrL7irn' n' (i) Pm' n' (j^) ~\~ ^ ^ Pm'n' (j^) 
rn'n' m'n' 



dt 



Sun'F,n,n'it)F^l{s)g';i'{t,s)C,t 



]{5mm'5i^m'Fl,n{t)Fi,l{s)g'^\ {t, s)Cst 

F^v{t)FUs)gl^{t,s)Cts) 



v 



(B22) 

From a computational point of view, this is the most 
useful form of the master equation. However, it is also 
useful to see the basis independent form of the master 
equation. Using the completeness relations of the basis 
states in the master equation above, one can remove the 
summations. After taking into account the possibility of 
multiple baths, one then ends up with the basis inde- 
pendent form of the master equation, which is given by 
Eq. (j66|) in the main text. 
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In order to solve the master equation, we first note 
that the system Hamiltonian is generally explicitly time- 
dependent. In such a case, we cannot, in general, cal- 
culate Us{t,s) in a simple way. Furthermore, solving 
the master equation itself becomes much more involved 
because of this explicit time-dependence, since the inte- 
grand in the master equation now depends explicitly on 
t. For our purpose, however, since we are only inter- 
ested in short times, we can still use a straightforward 
method to solve the master equation. Even though we 



know Us{t,s) for our case, we choose to start from the 
Hamiltonians - this serves as a check that the Hamiltoni- 
ans have been calculated correctly. We first use the split- 
operator method [63l465j to calculate Us{t,s), then per- 
form the integration in the master equation numerically, 
and we then finally solve the differential equation using 
the fourth-order Runge-Kutta (RK4) algorithm [gBI • We 
also note that there are considerably more involved and 
more efficient methods to solve such driven open system 
problems based on a decomposition of the bath spectral 
density (see, for example, Ref. [13] )• 
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